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ABSTRACT 

Propagation of torsional Alfven waves along magnetic flux tubes has been ex- 
tensively studied for many years but no conclusive results regarding the existence 
of a cutoff frequency for these waves have been obtained. The main purpose of 
this paper is to derive new wave equations that describe the propagation of linear 
torsional Alfven waves along thin and isothermal magnetic flux tubes, and use 
these wave equations to demonstrate that the torsional wave propagation is not 
affected by any cutoff frequency. It is also shown that this cutoff-free propagation 
is independent of different choices of the coordinate systems and wave variables 
adopted in the previous studies. A brief discussion of implications of this cutoff- 
free propagation of torsional tube waves on theories of wave heating of the solar 
and stellar atmospheres is also given. 

Subject headings: stars: atmospheres - MHD - stars: late-type - wave motions 

1. Introduction 

Direct measurements show that solar magnetic fields outside sunspots are concentrated 
into flux-tube structures located primarily at the boundaries of supergranules (e.g., Solanki 
1993). Similar magnetic structures are likely to exist in late-type stars (e.g., Saar 1996, 
1998). The fundamental modes supported by solar and stellar magnetic flux tubes can 
be classified as sausage, kink, torsional Alfven and fluting modes (for details, see excellent 
reviews by Hollweg 1990, Roberts 1991, and Roberts & Ulmschneider 1997, and references 
therein). The role played by these waves in the heating of different parts of the solar and 
stellar atmospheres was discussed by Narain & Ulmschneider (1996) and Ulmschneider & 
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Musielak (2003). The energy carried by sausage and kink tube waves was used as the input 
to the theoretical models of stellar chromospheres constructed by Cuntz et al. (1999) and 
Fawzy et al. (2002a, b). 

In many theoretical studies of tube waves, it is assumed that these waves are linear and 
they propagate along vertically oriented and thin flux tubes. Under these conditions, the 
waves do not interact with each other, so they can be investigated independently. Important 
studies of sausage tube waves were performed by Defouw (1976), who derived the wave 
equation for these waves and demonstrated that there is a cutoff frequency, which restricts 
the wave propagation to only those frequencies that are higher than the cutoff (see also 
Webb & Roberts 1979; Rae & Roberts 1982; Edwin & Roberts 1983; Musielak, Rosner, & 
Ulmschneider 1987). Similar studies of kink tube waves were performed by Spruit (1981, 
1982), who derived both the wave equation and the cutoff frequency for these waves (see 
also Musielak & Ulmschneider 2001). 

Propagation of torsional Alfven waves along solar and stellar magnetic flux tubes was 
extensively studied in the literature (e.g., Parker 1979; Priest 1982, 1990; Edwin & Roberts 
1983; Hollweg 1985; 1990; Poedts, Hermans, & Goossens 1985; Ferriz-Mas, Schiissler, & An- 
ton, 1989; Roberts 1991; Ferriz-Mas & Schiissler 1994; Roberts & Ulmschneider 1997; Hasan 
et al. 2003; Noble, Musielak, & Ulmschneider 2003). Two different approaches were consid- 
ered and different sets of wave variables were used. In the first approach, the propagation of 
the waves was described in a global coordinate system (e.g., Ferriz-Mas et al. 1989; Noble et 
al. 2003), while in the second approach a local coordinate system was used (Hollweg 1978, 
1981, 1992). The momentum and induction equations derived by Ferriz-Mas et al. (1989) 
were adopted by Ploner & Solanki (1999) in their studies of the influence of torsional tube 
waves on spectral lines formed in the solar atmosphere. In numerical studies of torsional 
tube waves performed by Kudoh & Shibata (1999) and Saito, Kudoh, & Shibata (2001), the 
basic equations originally derived by Hollweg were extended to more than one dimension 
and nonlinear terms were included. 

The specific problem of the existence or non-existence of a cutoff frequency for torsional 
Alfven waves propagating along thin and isothermal magnetic flux tubes has not been dis- 
cussed in the literature. An exception is the paper by Noble et al. (2003), who studied the 
generation rate of torsional tube waves in the solar convection zone and introduced the cutoff 
frequency, defined as the ratio of the Alfven velocity to four times the pressure (or density) 
scale height, for these waves. In this paper, we revisit the problem by deriving new wave 
equations that describe the propagation of torsional tube waves and demonstrating that this 
propagation is cutoff-free. We also show that the cutoff-free propagation is independent of 
different choices of wave variables and coordinate systems used by Ferriz-Mas et al. (1989) 
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and Hollweg (1978, 1981, 1992). 

The paper is organized as follows. The momentum and induction equations derived 
in a global coordinate system are presented in §2. In §3, new wave equations for torsional 
Alfven waves propagating along thin and isothermal magnetic flux tubes are obtained and 
it is shown that the wave propagation is not affected by any cutoff frequency. The fact 
that this cutoff-free propagation is independent of different choices of wave variables and 
coordinate systems is demonstrated in §4. A brief discussion of our results is given in §5, 
and conclusions are presented in §6. 



We consider an isolated and vertically oriented magnetic flux tube that is embedded 
in a magnetic field-free, compressible and isothermal medium. The tube has a circular 
cross-section and is in temperature equilibrium with the external medium. Let us introduce 
a global cylindrical coordinate system (r, 0, z) , with z being the tube axis, and describe 
the background medium inside the tube by the gas density po = po{r,z), the gas pressure 
Po = Po( r , z) and the magnetic field B = B or (r, z)f + B oz (r, z)z. The physical properties of 

— * 

the external medium are determined by p e = p e (r,z), p e = p e (r,z) and B e = 0. Moreover, 
we also have T = T e = const. 

To describe torsional Alfven waves, we introduce v = v^r, z,t)(f> and b = 6^(r, z,t)(f>, 
and assume that the waves are linear and purely incompressible, which means that both the 
perturbed density p and pressure p can be neglected. As a result of these assumptions, the 
propagation of the waves is fully described by the momentum and induction equations. The 
(/•-component of the momentum equation can be written in the following form: 



2. Basic equations 



dt\r ) 



47rp r 2 



1 



B 0r — + B 0z — (rbj) = , 



(1) 



and the 0-component of the induction equation becomes 





The derived momentum and induction equations are our basic equations for all the results 
derived and discussed in this paper. 
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3. Propagation along thin magnetic flux tubes 

3.1. The thin flux tube approximation 

Solar magnetic flux tubes are considered to be thin if their magnetic field is horizontally 
uniform, which means that at a given height all magnetic field lines have the same physical 
properties (e.g., Priest 1982). The essence of the so-called thin flux tube approximation 
(Roberts & Webb 1978, 1979; Spruit 1981, 1982; Priest 1982; Hollweg 1985; Ferriz-Mas, 
Schiissler, & Anton, 1989; Ferriz-Mas & Schiissler 1994; Musielak et al. 1995; Roberts & 
Ulmschneider 1997; Hasan et al. 2003) is that radial expansions around the axis of symmetry 
can be truncated at a low order. For the radial component of the magnetic field the leading 
term is of first order, 

~dB 0r (r,z)~ 



B 0r (r, z) = B 0r (r, z)\ r=0 + r 



dr 



+ , (3) 



r=0 



since Bo r (r = 0, z) — at the symmetry axis (cf. Ferriz Mas & Schiissler 1989). Away 
from the tube axis B Qr (r, z) can be expressed in terms of r and B 0z (z). Using the solenoidal 
condition V • B = 0, we obtain 



B 0r (r, z) = - r -B' 0z (z) , (4) 



where B' 0z = dB 0z /dz. 



The thin flux tube approximation also requires that p = po(z), p = po(z), B 0z = 
B 0z (z), p e = p e (z), p e = p e (z) and T = T e = const. In addition, the horizontal pressure 
balance must be satisfied, p + Bq z /8ti = p e at r = R t , where R t is the tube radius. The 
increase of R t with height is determined by the conservation of the magnetic flux ttR^Bq z = 



const. As a result of the above assumptions, the Alfven velocity ca = B 0z / a/47T/9 remains 
constant along the entire length of a thin and isothermal magnetic flux tube, and B' 0z = 
—B 0z /2H, where the pressure (density) scale height H is also constant. For reasons explained 
in the next subsection, the wave variables and b<f> are considered here to be functions of 
time and both spatial coordinates r and z. 



3.2. Wave equations 

Using Eq. (0J and taking into account the fact that the variables r and z are independent 
in the global coordinate system, we write Eqs. ([1]) and ([2]) as 

dv d B' 0z ( db^ \ B 0z db, 



at 8Tcp \ or J Anpo oz 
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and 



db <t> , B 'oz { r dv* 



dt 



dr 



-B, 



dz 



. 



(6) 



We combine the above equations and derive the wave equations for the wave variables 
v^r.z.t) and b^r, z,t) 



d 2 v 4> 2 d 2 v 4> 



dt 2 



dz 2 2H dz 1QH 2 



i r 



dv ri 



V<f, 



) (£) 



and 



d% 
dt 2 

(in) 



\AH) dr 2 



9% + 2 <9% 



1 dv d 



drdz AH dr 



= 



(7) 



-h dbj 



dz 2 2H dz 
d 2 bs . „ d 2 b, 



+ 2 



. 3 db d 



\AH J dr 2 drdz AH dr 



= 



(8) 



The derived wave equations show that the value of the parameter (r/AH) determines the 
contributions of the r-dependence of the wave variables and b<f> to the propagation of 
torsional tube waves. For wide flux tubes, the contributions are important, however, for very 
thin flux tubes with (r/AH) << 1, the contributions become negligible. It must be noted 
that the limit (r/AH) — * is not allowed because v</,(r, z, t)\ r= o = and b^(r, z, t)\ r= o = 0. 
From a physical point of view, this means that a flux tube reduced to a single magnetic field 
line cannot support torsional waves. 

Since the derived wave equations have different forms for and b^, the wave variables 
behave differently. In the following, we transform these wave equations to new variables that 
obey the same wave equations. 



3.3. Transformed wave equations 

Using the transformations v^r, z,t) = v(r^z^t)p~ x l^ and &^(r, z,t) = b(r, z,t)p 1 / 4 (see 
Musielak et al. 1995; Musielak & Ulmschneider 2001; Noble et al. 2003), we obtain 



dt 2 



d 2 v 
dz 2 



c A 



\AHJ 



d 2 v d 2 v 1 



dv 



(ah) (>r ~0r0z I// Or 



= 0, 



(9) 



and 



d 2 b 9 d 2 b 9 / r \ 

ci— — c • I — i 



dt 2 



l dz 2 



\ahJ 



r \ d 2 b ^ d 2 b 
AH/ W 2 + drdz 



J_db 
AHdr~ 



(10) 
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Clearly, the behavior of the transformed wave variables v and b is identical. 

To remove the first-order derivatives from the above equations, we use the transforma- 
tion d( = (AH/r)dr, which gives 

- 2c -aSc) wc - z ' t); 6(c - *< t)] = • 

where ( = AH In \r\. This is the most general equation that describes the propagation of 
torsional waves along thin and isothermal magnetic flux tubes. The equation shows that 
there is no cutoff frequency for torsional tube waves (see Sec. 3.4). 

3.4. Dispersion relation 

Since all the coefficients in Eq. (ITT1) are constant, we make Fourier transforms in time 
and space, and derive the following dispersion relation 

u 2 = (k 2 z + 2k z k c + kl)c\ , (12) 

where u is the wave frequency and k z and k^ are the z and £ components of the wave vector 
k, respectively. Note that the same dispersion relation is obtained for each wave variable. 

Let us define k = k z + k^ and write 

uj 2 = k 2 c 2 a , (13) 

which shows that the propagation of linear torsional Alfven waves along thin and isothermal 
magnetic tube waves is not affected by any cutoff frequency. 



4. Other approaches 

To demonstrate that the propagation of torsional tube waves is cutoff-free, we used the 
global coordinate system and the original wave variables and b^, which were transformed 
to the new variables v and b. Two different approaches were developed by Ferriz-Mas et al. 
(1989), who adopted the same coordinate system but used different wave variables, and by 
Hollweg (1978, 1981, 1992), who chose a local coordinate system and introduced different 
wave variables (see also Edwin & Roberts 1983 and Poedts et al. 1985). In addition, Noble 
et al. (2003) considered the global coordinate system and used the wave variables and b^. 
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However, their assumption that Bq t = was inconsistent with the solenoidal condition, which 
makes their claim of the existence of the cutoff frequency for torsional Alfveen waves invalid. 
In the following, we demonstrate that the momentum and induction equations derived by 
Ferriz-Mas et al. and Hollweg lead to the same results as those found in Sec. 3 of this paper. 



4.1. Global coordinate system and different wave variables 

In their work on propagation of waves along thin magnetic flux tubes, Ferriz-Mas et al. 
(1989) used the global coordinate system and derived the first and second-order equations 
that describe the propagation of sausage, kink and torsional Alfven tube waves. In their 
approach, each wave variable is expanded in a Taylor series and, specifically for torsional 
tube waves, the new variables v^i and b^i, which represent the first order expansion in the 
series, are introduced. These variables are given by 

v^{z,t) = ^| r=0 and b^(z,t) = ^-\ r =o ■ (14) 
Or or 

Since ity(r, z, t)\ r =o = and b^r, z,t)\ r=Q = (see Sec. 3.2), we may use Eq. (1T4"]) to 
write v$(r,z,t) = rv^z^t) and &^(r, z, t) = rb,pi(z,t) in first order. Substituting these new 
variables into Eqs. (jSJ) and ([6]), we obtain 

and 

db^ dv^ 

~W ~ B ^ = ° ' (16) 

which are the same equations as those obtained by Ferriz-Mas et al. (1989, see their Eqs. 
14 and 16; note that our B r0 corresponds to their rB rl , and v z and v r vanish in our case). 
Despite the fact that v$\ and b^\ depend solely on z and t, the above equations are not valid 
at the tube axis (see discussion above). 

The wave equations resulting from the above momentum and induction equations be- 
come 

^-4^0, (17, 

and 

d 2 hi od%i c 2 A db^ c\ 

~W ~ CA ~dz^ ~ H^z~ ~ 4fP b4>1 ~ ° ' (18) 
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Clearly, the derived wave equations have different forms, which implies that the wave vari- 
ables v^i and b^i behave differently. To remove the first-order derivative from Eq. ffT8]) . we 
use the transformation b^i(z,t) = b ( j,i(z,t)p 1 / 2 , and obtain 

^-4^ = 0. (19) 

Hence, the behavior of the wave variables and b^i is the same and there is no cutoff 
frequency that affects the wave propagation. This is an important result as it shows that 
the non-existence of a cutoff frequency for torsional tube waves is independent of the choice 
of the wave variables. 



4.2. Local coordinate system and different wave variables 

Propagation of torsional Alfveh waves along magnetic flux tubes can also be described 
by using a local orthogonal curvilinear coordinate system (£, 9, s), with s being the length 
measured along a magnetic field line, 9 the azimuthal angle about the axis of symmetry, and 
£ a coordinate in the direction £ = 9 x s. In this case, B = B 0s (s)s, B ^ = and B oe = 0. 
We also have v = vo(s,t)6, b = bg(s,t)9 and R = R(s), where R represents the distance 
from the magnetic field line to the tube axis. This approach was first considered by Hollweg 
(1978), who also applied it to solar magnetic flux tubes (see Hollweg 1981, 1992). 

Following Hollweg, Jackson, & Galloway (1982), the curvilinear scale factors are h$ = R 
and h s — 1, and we determine h% from the condition R B 0s = const, which results from 
V-I?o — 0. To conserve the magnetic flux, we must choose = R. Using these scale factors, 
the explicit form of the momentum and induction equations is 

d (v e \ B 0s d 

(Rb e ) = , (20) 



^(Rbe)-R 2 B 0s ^-( V 4)=0. (21) 



dt\RJ Anp R 2 ds 

and 

It is easy to see that the magnetic field B (s) in the local coordinate system can be 
described in the global coordinate system (see Eqs. [Hand [2]) by the B Qr (r, z) and B 0z (r, z) 
field components. This means that the following relation must hold between the spatial 
operators in these two coordinate systems 

d d d 

Bqttt + B 0z — = B Qs — . (22) 
or az os 
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This relation is consistent with the fact that Bo ■ V must be the same in the global (the LHS 
of Eq. [22]) and local (the RHS of Eq. [22]) coordinate systems. Moreover, v$ and vg are also 
related, as the former can be treated as a projection of the latter on the 0-axis of the global 
coordinate system. Obviously, the same is true for the wave variables and bg. 

We follow Hollweg (1978, 1981) and introduce the new variables x = vg/R and y = Rbg. 
The wave equations for these variables are 

w-Z& = °- (23) 

and 

^ 9 ( o^\ = , (24) 



dt 2 ds \ A ds, 

where in general ca = cu(s). However, for thin magnetic flux tubes ca = const (see Sec. 3.1) 
and the wave equations become 

|^-4^)KM);yM)] = o. (25) 

Again, no cutoff frequency exists. It is a significant (but expected) result that this non- 
existence of a cutoff frequency for torsional tube waves is independent of the choice of the 
cooordinate system and the wave variables. 



5. Discussion 

We considered the propagation of linear torsional Alfven waves along thin and isothermal 
magnetic flux tubes using the global coordinate system, and derived new wave equations 
describing this propagation. The derived wave equations were then used to demonstrate that 
no cutoff frequency exists for these waves, which means that torsional waves of any frequency 
are freely propagating along the tubes. We also showed that this result is independent of 
different choices of the coordinate systems and wave variables adopted by Ferriz-Mas et al. 
(1989) and Hollweg (1978, 1981, 1992). 

As first shown by Defouw (1976) for sausage tube waves and by Spruit (1981) for kink 
tube waves, the propagation of both modes is affected by their corresponding cutoff fre- 
quencies. With their cutoff-free propagation, torsional Alfven waves seem to be exceptional 
among the tube modes. In general, the existence of cutoff frequencies is caused by either 
gravity or gradients of the characteristic wave velocities, which result from an inhomogeneity 
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of the background medium. In the cases discussed in this paper, the characteristic wave ve- 
locities are constant for all tube modes because of the thin flux tube approximation. Hence, 
it is gravity which leads to the origin of the cutoff through either stratification (sausage tube 
waves) or buoyancy force (kink tube waves). 

The fact that stratification leads to a cutoff frequency for acoustic waves propagating 
in a stratified and isothermal medium was first demonstrated by Lamb (1908, 1911). Since 
sausage tube waves are essentially acoustic waves guided by the tube magnetic field, and 
since they propagate in a stratified and isothermal medium inside the tube, it is stratification 
of the background medium that is responsible for the existence of the cutoff frequency for 
these waves. 

The nature of kink tube waves is significantly different than sausage tube waves and 
yet it is again gravity that is responsible for the existence of the cutoff frequency for these 
waves. The main reason is that magnetic tension and buoyancy are the restoring forces for 
kink tube waves, and that the buoyancy force through its dependence on gravity leads to the 
cutoff frequency (e.g., Spruit 1982; Hollweg 1985), which is lower than that for longitudinal 
tube waves. 

Now, despite some similarities between kink and torsional Alfven tube waves, the main 
difference is that magnetic tension is the only restoring force for the latter. Since linear 
torsional tube waves have only purely axisymmetric twists in the ^-direction and show no 
pressure fluctuations, the twists are neither coupled to the gravitational force nor affected 
by stratification. As a result, no cutoff frequency can exist for linear torsional Alfven waves 
propagating along thin and isothermal magnetic flux tubes. 

The cutoff-free propagation of torsional tube waves may have important implications 
on theories of wave heating of the solar and stellar atmospheres. The theoretical models 
of stellar chromospheres constructed by Fawzy et al. (2002a,b) are based on the amount 
of energy carried by acoustic waves and by sausage and kink tube waves; these waves are 
generated by turbulent motions in the solar and stellar convection zones. The models point 
to a missing amount of heating for stars with high levels of activity. It is likely that the 
energy carried by torsional Alfven waves could be used, at least partially, to account for these 
"heating gaps" . Since there is no cutoff frequency for torsional tube waves, a broad spectrum 
of these waves is expected to be generated in the solar and stellar convection zones. The 
waves of different frequencies of this spectrum may transfer energy to different parts of the 
solar and stellar atmospheres. Hence, new studies are required to determine the efficiency of 
generation of torsional tube waves and their dissipation rates. 
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6. Conclusions 

We derived new wave equations that describe the propagation of linear torsional Alfven 
waves along thin and isothermal magnetic flux tubes, and used them to demonstrate that this 
propagation is cutoff-free. Our study also showed that the result is independent of different 
choices of the coordinate systems and wave variables used by Ferriz-Mas et al. (1989) and 
Hollweg (1978, 1981, 1992). 

Since the existence of cutoff frequencies for sausage and kink tube waves is caused by 
stratification and buoyancy force, respectively, our results clearly show that neither strati- 
fication nor buoyancy force affects the torsional wave propagation. The physical reason is 
that magnetic tension is the only restoring force for these waves. 

This lack of any cutoff frequency for torsional tube waves implies that a broad wave 
energy spectrum for these waves will be generated in the solar and stellar convection zones. 
The energy carried by the waves of different frequencies of this spectrum may be used to 
account for the "heating gaps" discovered by Fawzy at al. (2002a,b) and may also contribute 
to the wave heating of different parts of the solar and stellar atmospheres. 

This work was supported by NSF under grant ATM-0538278 (Z.E.M. and S.R.) and 
NASA under grant NAG8-1889 (Z.E.M. and S.R.). Z.E.M. also acknowledges the support 
of this work by the Alexander von Humboldt Foundation. 
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